Spin-stripe density varies linearly with hole content in single-layer Bi2201 cuprate 
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We have performed inelastic neutron scattering measurements on the single-layer cuprate 
Bi2+ccSr2-ccCu06+y (Bi2201) with x=0.2, 0.3, 0.4 and 0.5, a doping range that spans the spin- 
glass (SG) to superconducting (SC) phase boundary. The doping evolution of low energy spin 
fluctuations (< 11 meV) was found to be characterized by a change of incommensurate modulation 
wave vector from the tetragonal [110] to [100]/ [010] directions, while maintaining a linear relation 
between the incommensurability and the hole concentration, 5 ~ p. In the SC regime, the spectral 
weight is strongly suppressed below ~ 4 meV. Similarities and differences in the spin correlations 
between Bi2201 and the prototypical single-layer system La2-xSr x Cu04 are discussed. 



The relevance of charge and spin stripes to the phe- 
nomenology of hole-doped cuprate superconductors has 
been gaining currency in recent years. For example, in- 
triguing similarities in the transport properties of stripe- 
ordered cuprates and YBa2Cu306+ 2/ (YBCO) in a high 
magnetic field have been demonstrated [BE]- In fact, 
direct evidence for field-induced charge-stripe order in 
YBCO was recently obtained in a nuclear magnetic res- 
onance study [3 j. The impact of stripe order on the 
Fermi surface [4j [5] has been proposed as one possi- 
ble explanation for the appearance of quantum oscilla- 
tions [6j [7] (though there are challenges with the sim- 
plest approaches [8] ) . There is also evidence of closely re- 
lated nematic order [9] from Nernst-effect measurements 
on YBCO [10] as well as from spectroscopic imaging of 
Bi2Sr2CaCu208+ 2/ with scanning tunneling microscopy 
(STM) [IIJCE2]. 

Some of the most stimulating evidence comes 
from real-space imaging of electronic modula- 
tions by STM in Bi 2 Sr 2 CaCu208 +?/ [IMS] and in 
Bi2- x Pb x ST2-z^zCu06^- y [16]. These short-range 
correlations are found to have a period of approximately 
4a, where a = 3.8 A is the lattice spacing along a 
Cu-0 bond direction. The corresponding wave vector of 
the modulations, qsTM, is (|,0,0) in reciprocal lattice 
units (27r/a); gsTM is observed to decrease with doping, 
varying in the range of 0.3 to 0.15 [T5| fT6] . 

Identifying qsTM with the wave vector q co asso- 
ciated with charge stripe order in cuprates such as 
La2- a; Ba :E Cu04 [17] leads to a conundrum, as q co grows 
with doping (at least for hole concentrations p < 1/8 
[18]), opposite to the behavior of qsTM- When spin 
stripe order also occurs, antiferromagnetic spin corre- 
lations are modulated at q so = ^q co - It is often pos- 
sible to observe incommensurate (IC) spin fluctuations 



split about the antiferromagnetic wave vector Qaf by 
q<s ~ qso even when there is no significant stripe order, 
as in La 2 _ a ,Sr :c Cu04 (LSCO) [19] [20]. Among possible 
resolutions of the conundrum, it might be that the na- 
ture of stripe correlations is not universal among different 
cuprate families, or that qsTM measures something com- 
plementary to q co . 

In this Letter, we present the results of inelastic neu- 
tron scattering measurements of low-energy spin excita- 
tions in the system Bi2+^Sr2- x Cu06+ 2/ (Bi2201), demon- 
strating that 5 = \q s \ ~ P for 0.01 < p < 0.12. This be- 
havior is remarkably similar to that of I^-^Sr^CuO^ 
even including the rotation of q so by 45° for p < 0.06 
[19] [20]. These results provide strong circumstantial evi- 
dence that qsTM does not correspond to q co ; instead, it 
more likely corresponds to a nesting of antinodal states 
close to 2/cp, where kp is the nominal Fermi wave vector 
[T6] |2T] . This is not incompatible with a stripe origin, 
but would involve modulations along the charge stripes 
rather than perpendicular to them. 

While several variants of Bi2201 have been studied in 
the literature, we chose to work with Bi2+ x Sr2-^Cu06+ 2/ 
because it is possible both to vary the hole concentra- 
tion in a controlled fashion and to grow sufficiently large 
crystals with the floating-zone technique, as previously 
demonstrated by Luo et al. [22 . We prepared single 
crystals of Bi2201 with x=0.2, 0.3, 0.4, and 0.5. The 
actual concentrations of Bi and Sr were determined by 
inductively-coupled plasma (ICP) atomic emission spec- 
troscopy (AES), and the hole densities were determined 
by measurements of the Hall coefficient at 200 K, fol- 
lowing [23] ; the results are listed in Table 1. The cor- 
respondence between p and x is consistent with the pre- 
viously reported results based on measurements of the 
Fermi-surface volume by angle-resolved photoemission 



2 



TABLE I. Characterizations of the Bi2+c C Sr2-c C Cu06+y crys- 
tals. Elemental concentrations were determined by ICP-AES 
and hole concentration p was determined from Hall effect mea- 
surements. 



Bi 22 Sr lg Cu0 6+J 
7k0=llmeV, 7=70K 



Bi 24 Sr 16 Cu0 6+J 
fto)=4 meV, 7=40 K 



X 


Bi 


Sr 


Cu 


V 


0.2 


2.173(1) 


1.825(1) 


0.989(2) 


0.12(1) 


0.3 


2.282(3) 


1.717(2) 


0.992(4) 


0.09(1) 


0.4 


2.376(1) 


1.619(1) 


0.992(2) 


0.06(1) 


0.5 








0.01(1) 



spectroscopy [24 j. Based on magnetic susceptibility mea- 
surements, spin-glass-like behavior was observed below 
3 K for x = 0.4 [25] and below 4 K for x = 0.5; neither 
magnetic order nor diamagnetism were detected above 
2 K in the x = 0.3 and 0.2 samples. According to Luo 
et al. [22], the superconducting transition temperature, 
T c , is rsj 1 K at x = 0.2, rising up to a maximum of 9 K 
at x = 0.05, as shown in Fig. (3^b). The reduced T c in 
this system compared to La substitution for Sr is likely 
associated with structural disorder [26] . 

Most of the inelastic neutron scattering measure- 
ments were performed on thermal triple-axis spectrom- 
eter TOPAN installed at reactor JRR-3, Japan Atomic 
Energy Agency (JAEA). The typical collimator selections 
were 50 / -100 / -Sample-60 / -180 / , and the final energy was 
fixed at 14.7 meV. To reduce contamination from high- 
energy neutrons, a sapphire crystal was placed before the 
sample. A pyrolytic graphite filter was placed after the 
sample to eliminate higher-order neutrons. Additional 
measurements below 4 meV were performed at the cold 
neutron triple-axis spectrometers HER installed in the 
Guide Hall of JRR-3 and SPINS at the NIST Center 
for Neutron Research. For each composition, a couple 
of single crystals with total mass of 10-15 grams were 
coaligned and positioned so that the scattering plane cor- 
responds to (/i, /c, 0). Some results for the x = 0.4 crystal 
were reported previously [25] . 

For consistency, we will continue to index the scat- 
tering in terms of a tetragonal unit cell with a t = b t ~ 
3.81 A, although the symmetry is actually orthorhombic, 
with in-plane basis vectors along [110] and [110] corre- 
sponding to a Q and 6 G , respectively. Although we cannot 
resolve the very small orthorhombic strain, we can dis- 
tinguish the 6* direction by the presence of superlattice 
peaks (at ~ 0.26*) corresponding to the modulation of 
the BiO layers. We find that b* runs in a unique direc- 
tion in each crystal (i.e., there is little, if any, twinning), 
and we will see that this results in a unique orientation of 
the IC spin fluctuations in the more underdoped crystals. 

Inelastic neutron-scattering spectra for the x = 0.2 
sample (jp = 0.12) obtained for an excitation energy of 
huj = 11 meV and a temperature of T = 70 K are shown 
in Fig. [TJa)-(c). Scans A and B exhibit IC peaks split 
about Qaf in the [100] and [010] directions, while the 
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FIG. 1. Inelastic neutron scattering spectra in 

Bi2.2Sri.8Cu06+y x=0.2 at 11 meV, 70 K (a)-(c) and 
x=0.4 at 4 meV, 40 K (d)-(f). Insets show the IC peaks 
geometry and the scan trajectory. IC peaks are shown 
in (a) [100] and (b)[010] direction in x=0.2, and (d)[110] 
direction in x=0.4 sample, respectively. 



transverse scan C shows no structure. The pattern is 
identical to that observed in the superconducting phase 
of LSCO [27 and twinned YBCO [28] [29]; however, the 
intensity at this and lower energies is weak compared 
to that from LSCO for the same p and mass, measured 
under identical experimental setups. 

Related scans for the x = 0.4 sample (p = 0.06) are 
shown in Fig. [TJd)-(f); these were measured at Huj = 
4 meV and T = 40 K using the SPINS spectrometer 
with Ef = 5 meV. Here we see that IC peaks are only in 
scan A', which is along 6*, with no IC peaks along scan 
C, in the direction of a*. Similar scans at Huj = 1 meV 
are reported in [25], where the intensity is shown to fall 
off with temperature in a fashion consistent with mag- 
netic correlations. An earlier study demonstrated that 
the signal falls off in higher Brillouin zones, as expected 
for a magnetic form factor [30 . The observation of a lon- 
gitudinal IC splitting along a unique orthorhombic axis 
corresponds perfectly with the behavior previously found 
in the spin-glass phase of LSCO [20] [31] . From the unique 
orientation we infer that static order is likely, and it is 
strongly indicated by bulk susceptibility measurements 
to occur below 3 K [25]; however, we were not able to 
detect IC peaks in elastic scattering for any of the sam- 
ples. Of course, there is a substantial background in the 
elastic channel from nuclear diffuse scattering resulting 
from structural disorder, and that limits the sensitivity. 

To illustrate the variation of the spin correlations for 
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FIG. 2. Inelastic neutron scattering spectra at 4-6 meV and 
low temperatures in E^^Sri.gCuOe+y with x=(a) 0.5, (b) 
0.4, (c) 0.3 and (d) 0.2. Spectra were measured along [110] 
and [010] directions for (a,b) and (c,d), respectively. 



all four of our samples, Fig. [2] shows representative scans 
through Qaf obtained with thermal neutrons at excita- 
tion energies in the range of 4-6 meV. For x = 0.3, the 
orientation of the IC peaks is the same as for x = 0.2, but 
the splitting 5 is slightly smaller. For both x = 0.4 and 
0.5, the instrumental resolution is too broad to resolve 
IC peaks, but the narrower width for x = 0.5 suggests a 
smaller splitting. 

For quantitative analysis, we model the scattered in- 
tensity I(Q,u) ~ x"(Q,^)(l - e-^/^)- 1 with the 
formula, 

X "(Q, w) = x» exp [- ln(2)(Q - Qaf ± q*) 2 /^] . 

x"(w) and n correspond to the local spin susceptibility 
and the peak- width (half- width at half-maximum), re- 
spectively. In fitting the data, the model intensity was 
convolved with the instrumental resolution function and 
a linear background was included. For the x = 0.2 and 
0.3 samples, we defined q# to include (5, 0, 0) and (0, 5, 0). 
For x = 0.4 and 0.5, we set q s = (5/y/2, 5/y/2, 0); fitting 
the high-resolution data of Fig.[TJd) yielded S = 0.057(5). 
Since the fitted values of n for x = 0.4, 0.3, and 0.2 in 
Fig 1(d) and Fig 2(c,d) were comparable (~ 0.04 r.l.u.) , 
we fitted the x = 0.5 data with k fixed at 0.04 and with 
the assumption of IC peaks oriented as for x = 0.4. 
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FIG. 3. (Color online) (upper panel) Hole concentration de- 
pendence of the incommensurability 5 of low-energy spin fluc- 
tuations in Bi2+a;Sr2-xCu06+y (blue circles) compared with 
results for LSCO (gray triangles) QUEUES] and YBCO (open 
squares) [28, 29 Dashed line represents S=p. (lower panel) 
Electronic phase diagram of Bi2+ccSr2-ccCu06+y Spin-glass 
transition temperature is plotted by open circles. T c data 
(filled circles) are from 33 . Dashed lines are guide to the 
eye. 



The values of 5 obtained for Bi2201 from the fitting are 
plotted as a function of p in the top panel of Fig.[3j where 
they are compared with results for LSCO [27l|31][32] and 
YBCO |2E1 E21. We find that S w p in Bi2201, which ap- 
pears to be quantitatively identical to LSCO and quali- 
tatively similar to YBCO. Comparing with the spin-glass 
and nominal superconducting transitions indicated at the 
bottom of Fig. [3J the rotation of occurs between the 
spin-glass and superconducting phases, just as in LSCO 

The degree of similarity between Bi2201 and LSCO is 
a bit surprising, given that photoemission studies have 
indicated significant differences [34]. In particular, the 
chemical potential in LSCO remains rather constant with 
doping for < p < 0.12, while it shifts downward linearly 
with doping in Bi2201 [34]. There have been variety of 
models proposed to explain the doping dependence of 
S in cuprates [35l [36], not all of which involve stripes; 
nevertheless, if stripes are involved, then it is reasonable 
to expect that q co « 2q^, and hence |q co | « 2p for Bi2201 
with p < 0.12 based on the present results. 

The experimental frequency dependence of the local 
susceptibility x" (uo) is plotted in Fig. [4^ a). For the 
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FIG. 4. (Color online) (a) Energy dependence of the local 
spin susceptibility x"^) f° r Bi2+ccSr2-ccCu06+y with x=0.2, 
0.3, 0.4 and 0.5. Results at low temperatures below 13 K 
and high temperatures of 70 K for x=0.2 and 40K for x=0.3 
are plotted by closed and open circles, respectively, (b) Hole 
concentration dependence of x"^) a ^ 4 meV (triangles) and 
6 meV (squares). Dashed lines in (a) and broad line in (b) 
are guides to the eye. 



two samples near the spin-glass regime (p = 0.01 and 
0.06), the low-temperature x"^) * s large at low energy, 
consistent with proximity to an ordered state. For the 
more highly-doped samples (p = 0.09 and 0.12), x"^) 
is strongly reduced at low energy, and there is not much 
change when the temperature is raised somewhat (open 
symbols). The doping trend for the low-energy weight is 
illustrated in Fig. Qb). While it is qualitatively consis- 
tent with behavior observed in cuprates, the rapid loss of 
low-energy weight in this doping range is different from 
LSCO [37] . In particular, the results indicate that the 
system is far from static spin order of the type associate 
with the 1/8 anomaly in LSCO. 

Our result S ~ p in Bi2201 provides further evidence 
of universal behavior of spin correlations in the cuprates. 
The implication for possible coexisting charge modula- 
tions has a clear implication for the interpretation of 
STM studies. The doping dependence of qsTM is not 
consistent with the simplest stripe interpretation. The 
electronic modulations are more likely due to 2/cp-like 



modulations associated with the large antinodal density 
of states p~6j [21] . We are aware that qsTM does not pre- 
cisely match 2&f measured by photoemssion [38]; how- 
ever, this is not a problem as the STM modulations have 
maximum amplitude at bias voltages comparable to the 
pseudogap energy, so that the qsTM need not be deter- 
mined by states precisely at the Fermi level. It may be 
relevant that the antinodal pseudogap energy is also the 
scale on which antiferromagnetic spin fluctuations be- 
come strongly damped [39]. Hence, the low-energy spin 
fluctuations and the STM modulations appear to detect 
different electronic features with distinct energy scales. 
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